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ABSTRACT 


Ae preectical method for identifying linear time in- 
YVerliant sysitems om the basis of arbitrary input-output 
records is reviewed and extended to handle the case where 
tie sysitem is not wmitially im the zero state. The method 
is implemented using a digital computer program composed 
Of a mMumerical integration subroutine and a subroutine for 
solving overdetermined sets of linear alge basaaic <qacnenames 
severmal examples are presented to demonstrate the accuracy 


and present capabilities of the procedure. 
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I. INTRODUCTION 


A. THE IDENTIFICATION PROBLEM 

In order to apply any of the modern techniques of 
control system design one must first have a mathematical 
medeleoL the systemetosbe controlled. The form of this 
model will depend on the design methods to be employed as 
well as on the physical characteristics of the system. 
Since most of the theory on the analysis and design of 
control systems is based either on the state space or 
transform representation of systems the vast majority of 
mathematical models will consist of either a set of state 
equations or a transfer fume tuan . 

Once it has been dcecided what basic form the mathe- 
matical medel Sshewld take the problem of determining the 
numerical welues ef the parameters arises. Parameter 
values can’often be determined from the laws of physics 
and the data supplied by a manufacturer or obtained through 
testing. This is not always the case however. Occasionally 
the laws of physics become mathematically intractable or 
aCe ae VCimappliGcable. Qumute oftem the values of certain 
key parameters are not available. It is in these cases 
that the identification problem arises. 

Pacomion probliemingengineeringgis thateof determining 
the output of a system based on a knowledge of the system 
nmewel, input, sand initial conditions. The identification 


problem is similar to this but here the unknown quantity is 





the system model. The input and output are assumed to be 
known. For the purposes of this thesis the identification 
problem can be stated as follows: 

Given - a record of the input and output of a system 
Over some finite period of time, 

Find - a mathematical model and the numerical values 
of the model parameters in such a way that the 
model will accurately describe the behavior 
of the er | 

It should be kept in mind that the problem of identi- 
fying a system solely on the basis of input and output data 
(eiemso called “black box’” problem) is very rare in engi- 
neering. Even in the case where none of the model parame- 
ters are known one wiili more than likely have a fair idea 
of whether the system is linear or nonlinear, time varying 
Or time invariant, what the order cof magnitude of the domi- 
amt time Censtants is, and what types of inputs and out- 
Mitseane to be expected. For this reason most engineering 
identification problems fall into the "grey box" category. 
This distinction may seem trivial but virtually all iden- 
tification techniques rely heavily on knowledge of the 
characteristics and quantities mentioned above. 

Pretewon tie Comtrol systems literature on system 
identification is quite vast there are no Known identifica- 
tion schemes capable of handling all identification prob- 
lemissePrectively. Choosing a method suitable for a given 


problem can become a formidable task. A paper summarizing 


most of the common approaches to the problem of identifying 
lumped parameter systems has been published by Nieman, 
Fisher, and Seborg [1]. A good discussion of the indus- 
trial applications of various methods has been published 
Dye eykhotifi, et al. (i24« 

One common approach used in linear system identifica- 
ttom As sehat of obtaining the impulse response of the sys- 
tem. Mishkin and Haddad [3] have developed a technique 
fer finding the impulse response based on samples of the 
Syouen Guiput vamd 1ts derivatives. A technique for esti- 
mating the impulse response on the basis of noisy input and 
output samples has been developed by Levin [4] and Kerr and 
Seupen [S). Turin [6] and Lichtenberger [7,8] have used a 
meschedetilter to» obtain an identification. mae wee of wa 
Weave noise or binary test signal and crosscorrelation has 
been “suggested by Hak and McMurtry [9]. The nowse limita- 
tions of the sample approximation, matched filter, and 
Cresseeumetetion identifie@ationm techniques have been in- 
v¥esumpemed= by Lindenlaub and Cooper [10]. 

Another commoneapproach to the identification»preblem 
wemee determine the coefficients of the differential or 
difference equation which describes the system. Kumar and 
Srimmanr [11] have employed the method of quasilinearization 
with some success. -Nagumo and Noda [12] have developed a 
learning approach to the problem.” Bass [13] has developed 
a method which uses modulating functions and works well in 


the presence of noise. Astrom and Bohlin [14] have developed 


a statistically optimal method of determining differential 
equation parameters known as the "maximum likelihood method." 
A similar method which is not optimal but is considerably 
Simpler computationally has been developed by Peterka and 
Smek [15,16]. It is known as the “prior knowledge fitting 
method."' An algorithm for determining state variable models 
of sampled data systems has been proposed by Ho and Kalman 
[17]. The algorithm performs quite well in the presence 

of noise, and has been extended to continuously operating 
systems by Eldem [18]. 

Methods of identifying nonlinear and distributed param- 
etter systems ame usually lima@ted to specific types of sys- 
Pome or CO Specific types of nonlinearities. This is 
ungoubweedly due to the wide wariety of nenlinea@rities en- 
countered in physical systems and the difficulty of finding 
ee medel capable of Characterizing themyall. Shinbrot [19], 
Mowery [20], Fairman [21], and Bellman, pee: [22] have 
all approached the problem of identifying nonlinear systems 
by assuming a particular form of differential equation is 
capable of describing the system and then developing methods 
around the form of differential equation chosen. Another 
common approach to nonlinear system identification is that 
of representing a system by a suitable functional polyno- 
mil relating the inpwt and output. Hsieh [23] uses this 
approach and a steepest descent algorithm to solve the iden- 
tification problem. Similar approaches have been taken by 


sumeson (24), Bose [25,26], and Hubbell [27]. 


Identification methods vary widely with respect to how 
much must be known about the system before the method can 
be applied. Some identification techniques require that 
prior estimates of all system parameters be available. 

Many methods restrict the allowable system inputs to a 
family of testing functions such as steps or binary pulses. 
In general, the less that is known about cca System and 
the tighter the constraints on the kind of signals which 
mey be Gepliaecd as Gnputsethe more difficult it is to find 


emethed Gapable of accomplishing the identification. 


i. OBMECTIVES OF THIS PAPER 

iis paper will present a study of an identification 
technique Criginmally suggested by Diamessis [28]. It is 
designed to identify lumped linear time invariant systems 
but has been extended by Diamessis [29] and Wang [30] to 
tele Gertaie types Of nonlinearities. The technique re- 
quires a knowledge of the system input and output over some 
finite time interval as well as a rough estimate of the sys- 
tommoncer. the system input need not be restricted to a 
eilds> Of CteSting functions, it can be completely arbitrary. 

Wilke wseme rdentification techniques which require 
Eicmealculation of derivatives of the input and output, the 
technique to be presented requires only integrals of the 
input and output. The advantages of numerical integration 
over differentiation are well known. Since any zero mean 


ores Component on the input or the output tends to be 





attenuated preatly by the integration process the system 
identification can be more accurate than the raw data used 
to accomplish it. 

Tire weneanaes: of this thesis is divided into three 
impor sectrons. in Citapter II the theoretical developmént 
of the identification technique is given. A method for 
identifying the initial conditions of the unknown system 
icealse presented. Chapter III presents a method for im- 
plemenitrige tie techniques developed in Chapter I]. Par- 
Prculageattemelon 1S given to the choice of numerical 
methods and to efficient programming techniques. Several 
examples are presented to demonstrate the capabilities of 
the method. In Chapter IV several recommendations for 
further study are made. Conciusions concerning the accu- 
raey and present limitations of the technique under con- 
Sideration are also discussed. Following the conclusions 
a complete listing of all computer programs used in the 


thesis is given. 
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Pi. SIPeENTIPICATION BY MULTIPLE INTEGRATIONS — 


A. GENERAL APPROACH 

The development which follows is similar to the devel- 
opment given by Diamessis [28] in 1965. There are a few 
notable differences however. The development given by 
Diamessis is restricted to the case where all initial con- 
ditions are zero. This is a rather serious restriction 
Since it may be difficult or impossible to find a point 
where the system is in the zero state if the systems opera- 
ftom 1Senot to be disturbed. 7600 initial conditions will 
not be assumed in the development which follows. A method 
for solving for the unknown initial conditions will be 
presented. Diamessis proposed tne formulation of a uniquely 
determined set of linear algebraic equations with the model 
parameters as unknowns. This development will make use of 
overdetermined sets of linear algebraic equations with the 
model parameters and initial conditions as unknowns. The 
overdetermined set of equations will then be solved using 
the method of least squares. It will be shown that this 
results in amore accurate identification when the accuracy 
of the available data is limited and the order of the sys- 
tem 1s unknown. 

Any single input, single output, lumped parameter, 
linear, time-invariant system can be described by a linear 


ordinary differential equation with constant coefficients. 


JLo 


Mae basic form Of this equation 1S given in Equation (1) 


along with a set of initial conditions. 


n n-1 
d rit) + ane] d ae + eeeeete eevee + aoy (t) (1) 
dt dt 
m 
= 5b eh LUC. moe tec FH Ut) 
mS alee - 


wien 1 abtaal conditions ; 


Thea 3 ner 1 p 
d iG en yon i d u(0) y 1 


it 
Cc 


y (0) Ye u(0) 


where; 


u(t) = system input 


y(t) = system output 

Tiiemademeltication problem to be treated here consists 
of determining n, m, and the various coefficients of the 
differential equation on the basis of input and output 
records taken over some arbitrary time interval. The ini- 
tial conditions will be assumed to be unknown. 

Taking the Laplace transform of Equation (1) yields 


Equation (2). 


s"Y(s) + a, ys" UY(s) + creeeeeees + a ¥(s) (2) 


= bs U(s) eee + DEUS) 


net 
ee a a eeoenveeeeee ee 


RZ 





The g: Gocriicuemts accowmt for the contributiens ofthe 
mimewel conditions. Dividing Equation (2) by eee 


equivalent to integrating n+l times in the time domain. 











J Y S) + YG} + ef ee eevee My) 
S co ae: * 36 el (3) 
Us) ae UCs) 
Dn n-m+1 ° Do ont 
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a Be a5 + ccoeoeo eve eee ee + Bo ntl 


Taking the inverse transform of Equation (3) results in 


Equation (4). 


oe Ske eek 
[year+ ay f f y(t)dt? + (4) 
O O O 
tk ck 
eee + ay | eeentleece } y(tjyath*} 
O O 
vk vk 
= Db i *n-mtl---s f ieoydtee l= F ° 
O O 
ck bk 
ooo + by f eeentleee. [ aC dlredagagm 
O O 
mek ck ar 
5 Bn-] / i dt? + eee +- co | eentlee { atari 
O O O O 


Since the system is time invariant nothing has been lost by 


setting the lower limit on the integrals equal to zero. 


rs 





Rearranging terms in Equation (4) and placing all terms which 


g@ememd on u(t), (®t), or t a brackets results in Equation 


). 


(5) 
i t 
k k ty ty 
ate n+l y(tjat™*? | y (tides 
O O 
O O 
“k “k geal 
LDP euikel. u(t)dt tees 
O O 
“k oF *k *k 
Bese fy ened u(tjydt™ ™ + 4 g el. 
m O 
O O O O 
ty ot ty 
+eeet cay dt? = y(t)dt 
O O O 


Biecwrceemas Ol the Input and Output are assumed to 
be Known a linear algebraic equation with the system param- 
SieerS and imatial cemdition terms as unknowns can be formu- 
lated by performing the indicated multiple integrations 
from zero to some time ty: A set of 2ntm+l equations can 
be obtained by letting ty take on Zntmtl different values. 
Assuming that the equations are linearly independent it 
Pielenowspe possible to solve for the ntmt+tl differential 
equation coefficients and the n initial condition terms. 

So far nothing has been said of how n and m are de- 


Zmeencass  ineoretically it should be possible to use any 
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n' and m' greater than the actual order of the system under 
coudy. If the order of the system is n withm input coef- 


‘fPeients then» one would expect the following: 


a, = 0 fon O08 < 1 < n'-=n-1 } 
b. = 0 for ( 0 < 1 -< mtem=1 ) 
with n'>n and m'>m. 


The model should essentially reduce itself to the right order 
Be settingenonessential terms equal to zero. 

Unfortunately the situation is not quite this simple. 
Due to the finite accuracy of all experimental data the lin- 
oar comations will not haves an exact solution. Fom certain 
types of inputs the linear equations will not even be lin- 
cmiy mdependent. Tleese preblems can be overcomeste some 
extent by formulating more than 2n+m+l equations which are 
required. The overdetermined set can then be solved using 
Pnewiieene@d of least squames. If the limited accuracy of 
tie Experimental data can be attributed to round off errors 
ici iitteg@ratinp the data over a time interval much greater 
than the sampling interval should remove much of the uncer- 
toweeyean the linear equation coefficients. This is be- 
cause the roundoff process can be modeled as a zero mean 
white noise process. The integral of the noise will ap- 
proach zero as the period of integration increases. 

Even the measures mentioned above will not solve the 
problem completely however. Due to the finite precision 


Mieed to represent numbers in a computer and the iterative 


pls 





nature of the numerical methods used to solve overdeter- 
mained Sets of linear equations it 1s impossible to obtain 
zero as a solution for any parameter. If the parameter 
should be zero the numerical method will return a very 
small but nonzero value. Although this will result in an 
incorrect estimate of the system order the error will not 
Ben serious in most cases. This is because the small co- 
efficients of the terms which should nonexistent will make 
their effect negligible. Examples in Chapter III will 
demonstrate this point. 

Although the development in this section has assumed 
a differential equation model of the system the same results 
could have been obtained if a transfer function or state 
Variable mcdei had been assumed. A simplse rearrangoment cf 


terms in Equation (2) results in Equation (6). 


aes») bes. + ceoeeeeoe eee + by 


Ui Ss ) eT ceilcns 4 ajs + a, 





Sat In (6) 


By defining a few new terms, Equation (7) results. 





m m-l 
Df Gr) : ecs + c.4s reese oe (7) 
n n-1 
U(s) s + a_4s a S| a ay 
where K = b 
m 
Coe be / be wor Oli ams 


peesct OL State equations may be formulated in a similar 
fashion. One convenient state variable representation is 
fisenean cquatien (8). It is based on the phase variable 


form of system representation. 
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| 
l 0 | Se Xy 0 (8) 
) \ 
X5 0 0 \ X 4 0 
| \ 
a i ‘ , 
oe eee |, 
e 1 \ e 
0 ) 0 i 
ey ee Se ae Ss 
x a, ay coer a, xy K 
sh 
y = oe al e e e ® a 0 ® ® s ol rs 
. x 
n 


Note the simple correspondence between the terms in the trans- 


fer function and the terms in the state equations. 


B. IDENTIFYING INITIAL CONDITIONS 

Py ilanyeraemtitication problems it is désiraple to Gem 
panecemeagc system model with the actual system by exciting fhe 
System model with the same input data that was used in the 
Menertication. By comparing the output of the model with 
the output of the system a rough idea of the accuracy of the 
model can be obtained. This will not be possible however 
Um@hesoutic initial conditions at the beginning of the input- 
output record are all known. 

From Equation (5) it can be seen that when the linear 


algebraic equations are solved for the unknown model parameters 


Ly 








the g. initial condition terms are also found. By taking 
the Laplace Transform of Equation (1) and specifically 
writing in the contributions of the individual initial 
conditions a relationship between the g5 terms and the ini- 


tae conditions can be found. 


(o>) 
86 sal al - , an- 1 : Yo 
7 1 
8] 29 a / Yo 
7 
wf 
x 
a 
7 
7 
f n=2 
Bn -2 Do ge ir 
n-l 
i 
| Sn-1 te Jj Yo 5 
Cy C5 eee 1 0 ul U, 
GB I O 
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Unfortunately Equation (9) requires the knowledge of 
m-1 derivatives of the input function u. It may be neces- 
Samy to calculate m-1 derivatives of the input using numer- 
ical techniques. This could cause the model and system 
GCmeput to differ slightly at the beginning of the output 
mecord bt as thesnatural respense dies out the records 
should converge. It may be possible to avoid this diffi- 
eiay in Weanmy Cases by choosing the beginning of the input- 
emeaplt record at a point where the input is relatively 


constant. 


OF SEUMPLIFICATIONS WITH ZERO INITIAL CONDITIONS 

Pimieany pmeplems i1t°will be possible to exercise com- 
mere pcOltrol Over the input to the system under study. If 
fimmsme@essible to obtain an input-output record beginning 
WMimem the S¥YStem iS in the zero state it will be possible to 
meoity tne Laentification procedure. Since the g, terms 
Teer all be zero if the system is in the zero state they 
meeecd not be included in the formulation of the linear alge- 
pmate CqUuations. This will reduce the number of unknowns 
moneontm+ |! to nemtl. 

It should be noted that additional information can 
Cmaegene incorporated into the identification procedure in 
Guer tO Simplify the problem. For example if the steady 
State gain constant were known the number of unknowns could 
Demeecaweced by wone. It is usually a simple matter to tell 


MieGtetuaysyctem has poles or zeroes at the origin. This 
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piece of information can easily be used to simplify the 
identification procedure still further. As a general rule, 


the fewer the unknowns the more accurate the identification. 
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Til. IMPLEMENTATION 


Ae. NUMERICAL METHODS 

ie identification technique presented in Chapter II 
Game be brokén into two basic steps. The first step con- 
~iereceotepertormming the multipfe integrations of tHe input 
me output and forming the overdetermined set of linear 
eleebraic equations. The second step consists of solving 
the linear equations for the unknown model parameters. It. 
is a distinct advantage of the easement techmique 
Umer Study that each of these steps can be carried out by 
subprograms that are readily available in virtually all 
modern computer centers. 

Simepeonemeanebeswandled by any nuarerical imeegratieon 
subroutine. Although there are quite a few highly sophis- 
ticated numerical integration procedures available, trape- 
goidal integration will give better results in most 
applications. There are several reasons why this is true. 
fact Of all, most of the more complex integration tech- 
imees periomm poorly when the function being integrated 
PomcmscONntiiueus. Since control system inputs are often 
discontinuous and since such discontinuities are quite 
decanaple trem an identification standpoint, complex inte- 
gration techniques are usually undesirable. Even when the 
Functions to be integrated are continuous the slight im- 
provement in accuracy offered by more advanced methods is 
not enough to justify the tremendous increase in computation- 
al load associated with their use. 
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Step two, the solution of the set of overdetermined 
linear equations, is a classical problem in several fields 
of mathematics and engineering. Unfortunately most of the 
classical techniques for solving such problems are not 
Diactical. They tend to magnify the errors introduced by 
the finite precision of the Beiter tO the point where 
the solution is meaningless. Fortunately several modern 
methods are available that display more acceptable behavior. 
The method used in this paper was developed by Golub [31] 
in 1965. The basic approach is to triangularize the coef- 
ficient matrix by performing a Choleski decomposition. The 
decomposition is accomplished by applying repeated House- 
nemdenr transfommations [32]. Once the coefficient matrix 
has been triangularized the wnknowns can be obtained bv 
Baek substitution. The method is quite stable numerically 


aod 1S Capable of handling illconditioned coefficient ma- 


ze Gmce's. 


B. GCireAgteerslilCS OF GOOD INPUT-OUTPUT RECORDS 

The accuracy with which a system can be identified is 
Senomoly Cepencent On the input-output record used in the 
iiemtiticatvon. Since parameters are identified on the 
PewtemOn Meictr Citect on the OUtput it will be impossible 
EO -laciterty a parameter Piles ¢ its €@fect 1S meastneaple . 
I< thie Inpmwt to a system vas a frequency spectrum that is 
PoscvOrelccac tiletorm over the frequency range of interest 


Cee Tdemtinie¢auron will probably be very good. If the 


Cid 





frequency spectrum of the input is confined to a narrow 

band of frequencies the identification will probably be 

very pad. Tt is well kiiown that signals with sharp dis- 
continuities have a broader bandwidth than slowly varying 
continuous signals. For this reason input-output records dis- 
Paying discontinuities and rapid time variations should 

me chosen. 

If step or ramp imputs are used in the identification 
fie value of the initial conditions will have to be known 
Pacmimcorporated into the set of linear equations. Since 
mec initial conditions will usually be zero when these 
types of inputs are used this should not cause any diffi- 
menetes. ie reason for this difficulty lies in the nature 
tiem initiaiecencition coeriicient terms. THe inte~ral 
Ewetitcients OL these terms are steps, ramps, and higher 
order polynomials in time. If the input is a step or a 
mome the coefficients of several model parameters will also 
Bewsteps, ramps, and higher order polynomials in tiie. 
rete will therefore be a direct relationship between model 
Marameter coefficients and initial condition term coeffi- 
cients. This will result in the linear equations having an 
Misinite number of solutions due to the linear dependence 
between all the individual equations in the set. Step and 
ramp inputs must therefore be avoided when the system ini- 
tial conditions are unknown. 

Since analog system data will have to be converted to 


digital form a suitable sampling interval will have to be 
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chosen. Experimental results have shown that a sampling 
meme ten to one hundred times shorter than the shortest 
system time constant works quite well. Lower sampling 
names May mitroduGe imaccuracies in the location of high 


Peeduency polessand zeroes. 


c. PROGRAM DESCRIPTION 

iimonden ste cvaluave experimentally the characteristics 
of the identification procedure under study a set of digital 
Gomputer programs was developed. The main identification 
meopmam 15 a direct implementation of the procedure develop- 
edein Chapter II. Subroutine SYSTEM is a simulation program 
pect was Written tc generate input-output data for the main 
program to process. 

In the beginning of the identification program several 
important parameters are defined. NP and NZ are rough es- 
Pim@etes Of the number of poles and the number of zeroes in 
the system to be identified. KPTMAX is the number of san- 
miles peants available in the input-output record. Each sam- 
ple point will consist of the time (T), the input amplitude 
fy ane the output amplitude (C). IPTS is the number of 
Semple points that will separate successive linear equa- 
pilome. in Otter words, every time the total number of points 
Pao wines @ multiple of IPIS a linear equation will be gen- 
erated. The total number of linear equations that will be 
Beticwated 1S cqual to KPTMAX/IPTS. When all of the Linear 


equations have been formed subroutine DLLSQ is called. This 
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subroutine is an implementation of the Golub algorithm for 
fomvimameverdetermaumed sets of linear equations. 
subroutine DLLSQ returns the values of the system 
model parameters and initial condition parameters of Equa- 
tion (2). In order to find the poles and zeroes of the 
System the output of DLLSQ is fed into RTPLSB. RTPLSB is 
a polynomial root finder which .uses a combination of the 
Newton-Raphson and Seniceasc ors MetnodsstOwmmnd thespoles 
and zeroes of the system. 

The remainder of the main identification program is 
devoted to output. The results of the identification are 
given in both transfer function and state variable form. 
The state variable form is referenced to the format used 
wm kewation (8). 

SuUbreOmemicm or oleMunreads in the tramstfer functionjer a 
System and computes the system output based ona set of 
arbitrary initial conditions and an arbitrary input wave- 
form. Each time subroutine SYSTEM is called by the iden- 
tification program it returns three numbers, the time T, 
the input waveform amplitude R, and the output waveform 
amplitude C. Subroutine SYSTEM prints out the transfer 
function and state variable representation of the system 
Peder Simulating so that the accuracy of the identification 
can be determined. 

Input-output recrods obtained from physical systems 
ameeTanely accurate to more than three or four Signiticamt 


digits. The data generated by subroutine SYSTEM is therefore 


hie 





rounded off by subroutine ROUND before being passed to the 
identification program. The number of significant digits 

in the data returned to the identification program may be 

varied by changing the value of the parameter NA in the 


simulation subroutine. 


D. EXAMPLES 

The following examples demonstrate how the accuracy of 
an identification depends on the accuracy of the input- 
@uiepet record, the sampling period, and the input waveform. 
They also show how the order of the system can be deter- 
mined from a trial identification run using an estimated 
order greater than the actual order of the system. 

Eeemole one illustrates the relationship between the 
accuracy of the input-output record data and the resulting 
miemtrtication. Iineorder to Minimize the effect of other 
maecors all imitial conditions were set equal to zero, a 
step input was used, and n' and m' were set equal to n and 
mo oexemple ome demonstrates the fact that there is a di- 
rect (almost linear) relationship between the accuracy of 
input-output data and the accuracy of the identification. 
heme stmat even when the input-output record contained only 
CWwOwemonitficant digits the identification of the system 
poles. am zeroes was within about three percent of their 
Gxaet Values. 

Example two illustrates the relationship between the 


sampling period used with the input-output records and the 
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pie secimaecyeot the: identification, imput records contain- 
ing discontinuities or rapid time variations should be 
chosen. 

Examples four through ten demonstrate what happens 
mien nn’ awe@m' are preaver than’n and m. In each example 
an identification is performed using an n' and m' greater 
Phone and me Usine the informatizon obtained from this 
identification a new value of n' and m' is determined. 
These new estimates are then used to perform a second iden- 
Piel Cation. 

in examples fowr and five the error mm the estt#mate of 
m° was made larger than the error in the estimate of n'. As 
peuecult ef the relative values of these two errors’ all he 
emcess polles cancelled Wath @¢xeess zeroes’ But since thive 
Weme Miierme Excess zeroes than excess poles some excess zeroes 
Remeamned. “Nocée however that for frequencies lower than the 
Semplime Wate thie excess zeroes have little or no effect 
Smetnie™ Delavier of the identified system. Experiments have 
Sivewn Giiee excess zeroes that do notecaneel wrthwexcesis 
potes are always of a frequency comparable to or higher 
tiem tiemisemnpwing rate, Using this prancipite andeestrma- 
timpenew values for n' and m' results in an identification 
vien tas the correct number of poles and zeroes and is 
Momewaccurate than the first identification. 

in esempWes six, seven, and eight thewerror in the 
eseimate of n' was equal to the errer in the estimate of 


me eswa result, @1l excess poles and zeroes cancelled 
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with each other leaving a system of the correct order. 
Note that by repeating the identification with» the correct 
value of n' and m' it was possible to improve the accuracy 
Omithe identificatweon. 

In examples nine and ten the error in the estimate 
of n' was made greater than the error in the estimate of 
m'. As a result, all excess zeroes cancelled with excess 
poles but some excess poles remained. Note that for fre- 
Mienc ies lower than the sampling rate the excess poles 
have negligible effect. Experiments have shown that excess 
Bomlesethat do mot eamcel with excess zeroes are almost al- 
tiays Of a frequency Comparabie to or greater than the sam- 
meine rate. Using this primciple and estimating new values 
meen and i’ @resmited im a correct idemtification of ile 
Simulated systems. 

Experiments have shown that identifications involving 
Mie dalcecited xeess zeroes are generally more accurate than 
maentitications involving uncancelled excess poles. For 
eis reasom it 1s best to set.m' close to n' when identi- 
fying an unknown system. 

wee ene experimental findings described above a 
Simple procedure for determining n and m can be formulated. 
Mist, —80ess aan’ which 1s greater than the order of the 
system to be identified. This should not be too difficult 
Hemost cnpineecring identification problems. Let m' be 
equal to n' or m'-1. This will guarantee complete Games l= 


lation of all excess poles and zeroes or partial cancellation 
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Srecxcess poles and zeroes with excess zeroes remaining. 
After making a trial identification with the values of n' 
and m' mentioned above, estimate new values for n' and m' 
based on the reasoning in the examples. The new values 

of n' and m' should now be correct. By performing the 
identification with these values of n' and m' a good iden- 


freeacation should result. 
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EXAMEDE Ie 


STotemere BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
l -1 200000009000 0.0 
2 -10.99000099000 0.9 
3 -100.C000000000 0.0 

ZEROES REAL IMAGINARY 
1 -3.C0000009080 0.9 
2 —-30.00090CO00C0 0.0 

GAIN CONSTANT = 190.0000000CC0O 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 1600.0000000000 
2 1110.6009000050 
3 111.0200090900 
B VECTOR 
3 10.0009000000 
C VECTOR 
i 90.0000000000 
2 33 .C000000060 
3 1.6090C00009 


on 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 ~100.0901690225 0.0 
2 -0.9999970711 0.0 
3 -9.9999884241 0.0 

7 EROES REAL IMAGINARY 
1 ~2 49999937654 0.0 
2 -30.0000138488 0.0 

GAIN CONSTANT = 10.0000038147 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
i Ie FEMS U STIG 
2 1110.0003679024 
8 111.G6C01545177 
B VECTOR 
3 10.0000038147 
C VECTOR 
1 8929998545091 
2 33-0000076143 
3 1.0000C0000C0 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERICD 


0260600 SEC 
Ce20000 MSC 


it 


EQUATIONS FORMED = 129 
PRECISIGN OF DATA = Se DECTMAL PLACES 
INPUT FUNCTION = UNIT STEP 


or 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -~1090.0000291125 0.0 
2 -1.0000011522 0.0 
ZEROES REAL IMAGINARY 
1 -3 0000028353 0.0 
2 ~30.0000250480 0.0 
GAIN CONSTANT = 9.9999933243 
Syorc™ STADE VYAREPABIGES (PHASE FORM) 
A VECTOR 
1 1€00.6016C088045 
2 1110.0006141285 
3 111.8000319197 
B VECTOR 
3 9e 9999933243 
C VECTOR 
1 90.00C01602028 
2 33.0000278833 
3 1.Coccococcd 
PROGRAM PARAMETERS 
RECORD LEGNTH = 6.600890 SEC 
SAMPLING PERIOD = C.20000 MSC 
EQUATIONS FORMED = 120 
PRECISION OF DATA = ( DECIMAL PLACES 
INPUT FUNCTION . Ue T STEP 
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IOENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
l ~100.0006G 683416 0.0 
2 ~1.000011 8079 0.0 
LER@ES RE AL IMAGINARY 
1 —-3 20000272949 0.9 
GAIN CONSTANT = 9.9999847412 
SYSTEM STATE VARTABLES (PHASE FORM) 
A VECTOR 
1 10C00.0157822215 
2 1110.0653 743721 
3 111.9001130579 
BS Vecrur 
3 929999847412 
C MECTOR 
1 90.0010194754 
2 3320009941 70 
3 1.00CC0COC00 
PROGRAM PARAMETERS 
RECORD LEGNTH = 0.60000 SEC 
SAMPLING PERIOD = Oe 20000 MSC 
EQUATIONS FORMED = 120 
PECCISTGN GF DARA = 6 ‘BECTMAL ‘PACES 
INPUT FUNCTION = UNIT STEP 
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IDENTI 


FICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES 
1 95%, 
LE —te 
3 —Os 
EEROES 
1 26 
2 lr 


GAIN CONSTANT 


SYSTEM STATE VARIAB 


A VECTOR 
1 T9956 
Z 19% 
2 110. 
Bb VeerOR 
2 Ge 
C WECTOR 
1 89 « 
2 3% 
3 l. 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAPPLING PERIOD 
EQUATIONS FORMED 
PIRECISTON @F DAT 
INPUT FUNCTION 


REAL 

9968491528 
CO0001L7165 
9998390264 


IMAGINARY 
0.0 
0.0 
0.9 


REAL 
9999754422 
9992863605 


IMAGINARY 
0.0 
0.0 


946 S99S254 1 16 


CES (PHASE FORM) 


Seal S69 
94927167108 
9966898958 


9998254776 


OD OARS 1 
999 2 OMB 7 
O0COVVNOOGO 


0.60000 SEC 
0.20000 MSC 

120 

OS Cl wie eens leper > 
CNIS Sie 


i 


i 


A 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 —~99 69927835443 0.0 
2 -0.9991304951 0.0 
3 -9.9980 805636 0.0 
ZEROES REAL IMAGINARY 
l -~2-29981092096 0.0 
2 -29. 9981 688298 0.0 
GAIN CONSTANT = 9. 9954869232 
SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 
1 998. 8666303721 
2 1109 .6311321656 
cae 110. 9899946030 
B VECTOR 
3 9.9994869232 
C VECTOR 
1 89.9378461 845 
2 32 9962980388 
3 1-C000000000 
PROGRAM PARAMETERS 
RECORD LEGNTH = Ce6U00C SEC 
SAMPLING PERIOD = 0.2G0CC MSC 
EQUATIONS FORMED = 120 
PRECISION OF DATA = 4 DECIMAL PLACES 
INPUT FUNCTION = UNIT STEP 


56 


IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
l ~99.9937910594 0.0 
Z —~0. 9996622783 0.0 
ZEROES REAL IMAGINARY 
1 ~2 9996597024 0.0 
2 —~30. 0061026888 0.0 
GAIN CONSTANT = 9.9986057281 
SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 
1 999. 7584771470 
2 1110.0544578539 
3 110.9950 366432 
B VECTOR 
3 9.9986057281 
C VECTOR 
l 90 .0678270058 
2 33.0057533912 
3 120000060960 
PROGRAM PARAMETERS © 
RECORD LEGNTH = 0.60000 SEC 
SAMPLING PERIOD = 0.20000 MSC 
EQUATIONS FORMED = 120 
PRECISION OF DATA = 3 DECIMAL PLACES 
INPUT FUNCTION = UNIT STEP 


Si 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


BOLES REAL IMAGINARY 
i —-96'. SC1O45 76523 0.0 
2 -0.9763637189 0.0 
3 -9.8862584052 0.0 

ZEROES RE AL IMAGINARY 
1 —2 29406610975 0.0 

GAIN CONSTANT = 929451894760 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 953.6797208413 
2 Wee. 8846223636 
3 169. 6630799769 
B VECTOR 
3 929451894760 
C VECTOR 
1 871897913926 
2 32.5503855989 
3 1.0000000000 


PRCGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 


0e-60000 SEC 
0.20060 MSC 


EQUATIONS FORMED = ho 
PRECISION OF DATA = 2 DEC THA PIRACES 
INPUT FUNCTION = UNIT See 
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EXAMPLE 2e 


Sy ShaMeTOeBE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POWES 
1 
2 
5 


ZEROES 
1 
2 


GAIN CONSTANT = 


REAL IMAGINARY 
-1.0000000000 0.0 
-10.0000000000 0.0 
~100.0000009000 0.0 
REAL IMAGINARY 
-—3.-000C0COCGOCD 0.0 
~30.00000000C0 0.0 
10. 0009000000 


Siores Ssmave VARIABLES (PHASE FORM) 


A VECTOR 
1 
2 
3 


B VECTOR 
3 


COVECTOR 
i 
2 
3 


1000. 0000S00GC0 
11190.COCO00000C 
111.0900000009 


10-C00000009C 


90.COGN699000 
33 -0000CO00CO 
1. 0000000900 


a 





IDENT I 


SYSTEM TRANSFER FUN 
POLES 

1 ~101. 

2 ~le 

3 -~10. 
LBRGES 

Me =3'. 

2 -~30. 


GAIN CONSTANT = 


SYSTEM STATE VARIAB 


A VECTOR 
1 1010. 
2 iieeoliee 
3 ll2. 
B VECTOR 
& 10. 
C VECTOR 
1 90 
2 33. 
3 le 


PROGRAM PARAMETERS 
RECORD LeGNTH 
SAMPLING PERTOD 
EQUATIONS FORMED 
PRECISION OF DAT 
INPUT FUNCTION 


FICATION OF UNKNOWN SYSTEM 


CTION 

REAL IMAGINARY 
0610717819 0.0 
0000331312 0.9 
0017011716 6.0 

REAL IMAGINARY 
0002183624 0.0 
0553440036 0.0 

10.0887174606 


LES (PHASE FORM) 


8161284948 
8490926385 
0628060847 


0887174606 


1725949663 
0555623654 
0000000000 


0.66000 SEC 
0.19000 MSC 

120 

3° DECIMAL "PLACES 
UNTT Ser 


A 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1. -105.3244028971 0.0 
Ze -1.0000069733 0.0 
3 -10.0258525504 0.90 

LERDE S REAL IMAGINARY 
1 ~3 20605906385 0.0 
2 —30 2 4fehO561 855 0.0 

GAIN CONSTANT = 10.3944911957 


Svoueh STARE VARTABLES (PHRASE FORM) 


Ah WEC TOR 
1 5 5s SitGe2 9613 1.92 
Z TTT lap oo Po o2 211 
3 116. 3502624208 
B VECTOR 
3 10.3944911957 
C MECTGR 
1 Mie 43 h659562 
2 3344716468241 
3 1.0009909C90C0 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 


i 


1.2@@C0 SEC 
G.200G0 MSC 


EQUATIONS FORMED = 2.0) 
PRECISION OF DATA = 3 DECIMAL PLACES 
DNPUT FUNGiIOREN = UNIT STEP 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


P-OWES 
] —10r. 
2 -l. 
3 -10. 
LERGES 
] =ai6 
2 -30. 


GAIN CONSTANT = 


REAL 
Gate clos9 
0000162038 
60077785631 


IMAGINARY 
0.0 
0.0 
0.0 


REAL 
0001692909 
3047449110 


IMAGINARY 
0.0 
0.0 


FOr COS 13306 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
i 1C77. 
2 1194. 
g 118. 
B VECTOR 
3 10. 
C VECTOR 
] 90. 
Z S18) 
3 le 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 


EQUATIONS FORMED 
PRECIST@N OF DATA 


INPUT FUNCTION 


2489572831 
8806091150 
6472169649 


6655513306 


epee tess ©71LY 
3049142009 
C0C0900000 


I 


3e000G0 SEC 
0.50660 MSC 

120 

3 DEC PMAL PLACES 
UNIT STEP 
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EXAMPLE 36 


Srovene1 Oo GCE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0000090000 0.0 
Zz -10.0C00000000 C.0 
3 -~100.00000090000 0.0 

ZEROES REAL IMAGINARY 
1 -3 20000090000 0.9 
Z -30.000000090C 0.0 

GAIN CONSTANT = 19.0000600900 


SyourRM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 1000.C90C0906C0 
Z 1110.0000000V0C 
3 111.0900000000U 
beveGrer 
3 10.000000CO00GOG 
C VEGTOR 
1 90.-0900000000 
2 33 -QCO000N090C 
6) 1.00000900C0 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 ~99.6074 714180 0.0 J 
2 ~1.6627466213 0.0 
3 —10.0526991957 0.0 J 
ZERGES REAL IMAGINARY 
1 -3.110°9954642 0.0 
GAIN CONSTANT = 9. 9883327484 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 1064. 15364%23652 
2 1117.8649236164 
3 110.7229172350 
B VECTOR 
3 99883327484 
C VECTOR 
1 93 .0004760964 
2 330051168692 


3 1.0000006000 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 
EQUATIONS FORMED 
PREGLSTONe@F DATA 
INPUT FUNCTION 


6.60000 SEC 

0.10000 MSC 

120 

3 DECIMAL PLACES 
GAUSSIAN- SIGMA=0.010 SEC 


i 
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IDENTIFICATION OF UNKNOWN SYSTEM 


a 


Svs Tem TRANSFER FUNGTION 


P.GiIRES REAL IMAGINARY 
l ~98.99033657608 0.0 J 
2 ~0.9389787239 0.0 
LERGES REAL IMAGINARY 
] ~2¢ 8677895731 0.0 
2 ~29.6513565567 0.0 
GAIN CONSTANT = 949 DOB 285522 


Soren SiAME VARTABLES (PHASE FORM) 


A VECTOR 
ul Gieesea8 27259 
2 1974.4957479197 
3 109.6741207397 
B VECTOR 
3 eae 285522 
C VECHOR 
1 SSe033 8511612 
2 3225191461297 
o 1.0000090900 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 
EQUATIONS FORMED 
PmeCISIGN OF DATA 
INPUT FUNCTION 


0.60600 SEC 

Oc 1006 MSC 

120 

2.DECTRMAL, "PEACES 
GAUSSITAN- SIGMA=0-030 SEC 


tt 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 

1 -99.1859265460 0.0 o 
ZEROES REAL IMAGINARY 

l =2'> 18659890 38 C.0 

2 -29.64122032416 0.0 
GAIN CONSTANT = 9.9813451767 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 TIP SSCTOCo tors 
2 1026.2057811417 
3 109. 4592852321 
B VECTOR 
2 929813451767 
C VECTOR 
: 1420793075760 
Z 3129308621454 
3 1.000G0CCOCO 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 
EQUATIONS FORMED 
PRECISION OF DATA 
INPUT FUNCTION 


0-60C00 SEC 

0O.100CG MSC 

120 

a DEGUMAL PEAGES 
GAUSSIAN- SIGMA=0.050 SEC 
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EXAMPLE 4 
SyYoreMe TG BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 ~5.0600000000 

2 -~20. 00G00C00C0 
ZEROES REAL 
GAIN CONSTANT = 300.-0CGCOOG900 


SyoTEMMSTATE VARIABLES (PHASE FORM) 


A VECTOR 

1 L100-90CO000COO0N 

2 25-.G0S009C0C0 
B WECTOR 

2 309. OOICINOCOOO 
C VECTOR 

1 1. Q0000C0GCC 


INITIAL STATE VECTOR 
i 1.-.C0009000000 
Zz 0.0 


47 


IMAGINARY 
0.9 
0.0 


IMAGINARY 





IOENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL 
1 ~20. 1020603504 
2 (7064280895557 
3 ~0.4280895557 
4 1.90489471395 
> -4.9974142741 
ZEROES REAL 
Ls 21168. 5357639489 
2 -0.4681078779 
8 -0-. 4681978779 
4 1.0470216652 
GAIN CONSTANT = -0.9141926892 


Seeouern STADE VARIABLES (PHRASE FORM) 


A VECTOR 


mn & W %M 


B VECTOR 
5 


C VECTOR 
i 


mn f-& WwW NH 


-4826. 300213 
33052920433 
Prue. Oon) OF 
140.5228441 
24290670660 


—O% l49268921E-C1 


19201812159 
=o 256062005 / 
2399-647541 
—Z1166. 64657 
1.0c0090000 


48 


IMAGINARY 
C.0 

—-6- 1540893481 

627540893481 
0.0 
0.0 


IMAGINARY 
0.6 

=6.6 1 €82935 708 

6. 7682935708 
0.0 


Qj. a G&G G G& 


a Qa Qa Ge 





IDENTIFICATION OF UNKNOWN SYSTEM 


INITIAL CONDITION (G) VECTOR 


1 O. $983029652 
2 25eOUNSse179 
3 37.55868385 
4 TNS 123299 
5 =Wo2seo39G00 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 
EQUATIONS FORMED 
PRECISION OF DATA 
INPUT FUNCTION 


3.00000 SEC 
0.50009 MSC 

120 

2, OECI MAL PLACES 
PIECEWISE CONSTANT 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 —-4.9994757762 

2 ~20.0118644519 
ZEROES REAL 
GAIN CONSTANT = 300.1457519531 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 100.0488316 

2 25201134023 
8 WEGTIOR 

2 36021457520 
C VECTOR 

1 1.600000090C 


INITIAL CONDITION (G) VECTOR 
1 12001249200 
2 ZSZCOUSS0377 


50 


IMAGINARY 
0.0 
0.0 


IMAGINARY 





EXAMPLE 2S 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 -2.00900CC000 

2 -2.00000CC090 

3 -~20.00009000G0 
ZEROES REAL 

1 ~8.000000000C 
GAIN CONSTANT = 160-060C000000 


SYouEwe STATE WARTABLES (PHASE FORM) 


A VECTOR 
1 260.00000C0000 
2 93.00009000C0 
3 24.000000C000C 
B VECTOR 
2 160.cC000000000 
C VECTOR 
l 8. OCO0000C00G 
2 1.-0000000000, - 


INITIAL STATE VECTOR 


Ls 2-00000000C0 
2 1-0000000000 
3 0.0 


sph 


IMAGINARY 
320000000000 J 
-~3.COCOO00G000 J 


0.0 J 
IMAGINARY 
0.0 J 





SYSTEM TRANSFER FUNCTION 


POLES 


wm & W NM 


ZERGE S 
1 


@ 
5 
4 


GAIN CONSTANT = 


Seomet STATE VARIABLES (PHASE FORM) 


A VECTOR 


1 


wm ff W DM 


B VECTOR 


5 


C VECTOR 


1 


wm & W NY 


REAL 


i SI eee n2 Ie ho a Or | 
120012938647 
1.0012938647 

~1.9998854338 
~1.9998854338 


REAL 


~7-6D389C2797 
447.1906330529 
be29222266906 
lwee 92 22 26696 


2UBG3 6 G7 13 
69964430945 
1994.127061 
124. 3802348 
22025117861 


—~9- 3623618484 


~283455. 3928 
~27855. 70425 
ellos 740091 
~442.1711881 

1.000009000 


A 


~0. 3623618484 


IDENTIFICATION OF UNKNOWN SYSTEM 


IMAGINARY 
0.0 
-8.28128722093 
8- 8128722093 
—Je0G1 | Oooo Jo 
3. 0017339895 


IMAGINARY 
0.0 
0.0 
-9.0382447685 
920382447685 


aj. c& Ga Ga G 


a. GQ a Ga 





IDEN 


INITIAL CONDITION 
1 


wm ££ W NHN 


PROGRAM PARAMETER 
RECORD LEGNTH 
SAMPLING PERIO 
EQUATIONS FORM 
PRECISION OFsD 
INPUT FUNCTION 


TIFICATION OF UNKNOWN SYSTEM 


(G) VECTOR 
16.99753824 
388.4269451 
1684.166654 
30704.56626 
95381. 34714 


S 


3200000 SEC 
0.50000 MSC 

ZU 

2 DEG TMAL PLACES 
PIECEWISE CONSTANT 


D 
ED 
ATA 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 eel 5161 423 73 

2 -1.999/989237 

3 —108999 1989237 
ZEROES REAL 

1 ~8.03356/7123 
GAIN CONSTANT = 16024525451 660 


Sroven STATE VARTABLES (PHASE FORM) 


A VECTOR 
a 2615830183 
Z 93252427869 
S 24213121208 
B VECTOR 
3 160.4525452 
C WECTOR 
1 82033567712 
2 1.09000C00900 


INITIAL CONDITION (G) VECTOR 


1 Oe 99920347 
2 cece oc | 16 
] 1168.3 76804 


394 


IMAGINARY 
6.9 J 
—-3.0007T 298212 J 
3e0011201212 J 


IMAGINARY 
0.0 J 





EXAMPLE 6 


Syoueey TOMBE ITOENTIFIED 


SYSTEM TRANSFER FUNCT I ON 


POte S REAL IMAGINARY 
1 -3.9000000000 0.0 
2 -45.-000000CCO0S 0.0 

ZERGES REAL IMAGINARY 
1 -15.90000000CG0 6.0 

GAIN CONSTANT = 10.00G000G00C 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 135.0000000099 
2 48.0000000090 
B VECTOR 
2 1¢.0000990000 
Ce veeuaR 
1 15.00009000000 
2 1.000C000000 


INITIAL SFATE VECTOR 
1 1.0090000000 
2 0.0 


a0 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES 
1 


mn ££ WwW N 


LPERGES 
i 


2 
5 
4 


GAIN CONSTANT 


REAL 
~44,9825 156448 
Ue 8248831729 
0.8248831729 
=2e 97S 94uIao D 
~4.9057021934 


REAL 
Os 8 3629742371 
0.8362974237 
ioe 0 59007 3603 
~4.901 7160580 


= pe 7og2181933 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 


wn fF W NH 


B VECTOR 
5 


C VECTOR 
1 


wm & W WN 


OOM 2 «59S 9 
104719.1699 
LO ire TC ODd 
568.8244217 
51.24110636 


De Peg2(S7193 


21041.07999 
woes J05 59 
B20 749073 
13026812859 
1.0C00C009C 


56 


IMAGINARY 
0.0 
-16.880278213C 
16. 88027821 36 
0.0 
0.v 


IMAGINARY 

—16. 813963560) 

16.8/739835601 
0.9 
0.0 


qe. G&G. Ga GQ G@ 





IDENTIFICATION OF UNKNCWN SYSTEM 


INITIAL CONDITICN 


wm + W NH 


PROGRAM PARAMETER 
RECORD LEGNTH 
SAMPLING PERIO 
EQUATIONS FORM 
RRECTSTON OF BO 
INPUT FUNCTION 


(G) VECTOR 
14.99958219 
632200.0035658 
6068.388958 
183314.0231 
Ss She sre iS Sie! 


S 


Meoeoo> SEC 
Wa2e22Z ASC 

120 

4°DECIMAL PLATES 
PIECEWISE CONSTANT 


D 
ED 
ATA 


37 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSEEROEUNCT GON 


POLES REAL IMAGINARY 
2 -44.9468936931 0.90 

ZEROES REAL IMAGINARY 
1 -14.9948912787 0.0 

GAIN CONSTANT = 929912204742 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 134.8316990 
2 47.94669385 
B VECTOR 
2 929912204 74 
C VECTOR 
1 14.99480128 
2 1 -000900000 


INITIAL CONDITION (G) VECTOR 
1 W45998 46932 
2 584.3008098 _ 
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EXAMPLE qT 


SYSteeeaGaee IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


PUCES REAL 

1 -1.000000G000 

2 -1.-0000000009 

3 ~59.0C000CU000 
ZEROES REAL 

1 -~5.0CCGO0COCO 

2 ~20. 00CCO00CC0G 
GAIN CONSTANT = 5-000009000C000 


See eo tAtE VARIABLES (PHASE FORM) 


A VECTOR 

1 100.0099009006 

2 102.0099900000 

3 52.0C00CGCOUG 
B VECTOR 

3 5.0000000060 
C VECTOR 

} 100.,00000C00K - 

2 25. 0060000000 

3 1. 00C00G0C00 


INITIAL STATE VECTOR 


JL 2-00000G00G60 
2 1-909000C0CO 
3 G.0 


59 


IMAGINARY 

-1.c0occoscoo0cg 

1. 00000C0GCO 
C.0 


IMAGINARY 
0.0 
0.0 


J 
J 
J 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


P@LES REAL IMAGINARY 
1 me 9911920718 0.9 
2 121558667717 =236 leu ooo 
3 ie t 5506o7 7 17 2324137133914 
4 —C.99972239C9 —G997465371'5 
5 ~06:9997223909 0.9974653/18 

ZEROES REAL IMAGINARY 
1 1.21/9838023 “2305 (oo lace 
2 1-2179838023 23031 Seo lead 2 
3 —-19. 8075068598 0.9 
4 SS eehoZe PS012 C20 

GAIN CONSTANT = 429560489655 


Semen Stale VARIABLES (PHASE FORM) 


A VECTOR 
1 54739. 78240 
2 55794. 37204 
3 DBAS 4Ne5333 
4 53162985138 
5 49. 67890331 
B VECTOR 
5 40956048965 
C VECTOR 
1 55512. 60354 
2 e708. 114537 
3 58845004083 
4 22. 48682083 
5 16060090000 


60 


qj a. Ga G&G G& 


&. Gee. Cc 





IDEN 


INITIAL CONDITION 
1 


uw Ff W NY 


PROGRAM PARAMETER 
RECORD LEGNTH 
SAMPLING PERTIO 
EQUATIONS FORM 
Peels RON OF OD 
INPUT FUNCTICN 


TIFICATION OF UNKNOWN SYSTEM 


€G) VECTOR 
225.00612 73 
10975.84397 
117558. 8284 
6270060.545 
PiZs 2 c6.> 1 


S 


1220000 SEC 
O.2C00C MSC 

129 

4 DECIMAL PLACES 
PIECEWISE CONSTANT 


D 
ED 
ATA 
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IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -49.9862042435 0.0 J 
2 -0.9995303000 -1.00C 7669686 J 
3 -0.9995303000 1.0007669686 J 

ZERGES REAL IMAGINARY 
1 -4.9622611725 0.0 
2 -20.2671062457 0.0 

GAIN CONSTANT = 4.9685173035 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
sy 100.0021676 
2 101.9269468 
3 51298526484 
be VECTOR 
3 42968517303 
( VECTOR 
1 160.5736518 
Z Pree eCTIONHE. 
3 1.-000C00G96 


INITIAL CONDITION (G) VECTOR 


1 22520065829 
2 11494. 69251 
3 20485. 93332 


62 





EXAMPLE 8 


Syoner) nee oe IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


PGGES REAL 

] ~3.90000000090 

2 -8.0000000000 

3 ~40.CO000G0000 

4 ~40. GO0G000CCO 
ZERGES REAL 

1 -20.000GG60cCC00 
GAIN CONSTANT = 169.-00000G60000 


Peotem STATE VARTABLES (PHASE FORM} 


A VECTOR 
1 39C60.00200C00G0 
2 19795, 0C000C00GU 
3 2529. 00000C000C 
4 91.00000C0060 
B VECTOR 
4 160.0000000000 
C VECTOR . 
1 20.00000000C0 
2 1.00000000U0 


INITIAL STATE VECTOR 


1] 3-0000C00000 
2 220000000000 
2 1.00000C000G 
4 0.0 


63 


IMAGINARY 
0.0 
0.0 

~5. 9000006006 

5-« OOCO0COCLO 


IMAGINARY 
0.0 





IDENTIFICATION OF UNKNQWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES 
1 


wm & W NM 


ZEROES 
l 
2 


GAIN CONSTANT 


SYSTEM STATE VARIABLES 


A VECTOR 
1 


nm fF W NH 


B VECTOR 
5 


e VECTOR 
i 
2 
3 


REAL 
~40.0576848888 
-40.0576848888 

~8.0002609098 
=1.4530 33153 
—~22999 8483817 


REAL 
-1.4492 282436 
-—20m5 685168935 


56181. 06719 
67897.28681 
eos 35114 
2665.954710 
92. 56849238 


157.7509613 


29% O1E7T1692 
21.81789514 
1.G90°00000 


64 


57. TSOFECWS03 7 


Ve eASe FORT) 


IMAGINARY 

-4. 8648939727 

4 8648939727 
0.0 
0.0 
G.0 


IMAGINARY 
0.0 
Oe 0 


a a. Ga a Ga 


IDENTIFICATION OF UNKNOWN SYSTEM 


INITIAL CONDITION 


wm & W NY 


(G) VECTOR 
62.00108720 
5780206688 
169308.2337 
14097642537 
1708815.026 


PROGRAM PARAMETERS 


RECORD LEGNTH 
SAMPLING PERIO 
EQUATIONS FORM 
PRECISIUN OF D 
INPUT FUNCTION 


1648842 SEC 
Ce.24807 MSC 

120 

See CUMALOPIEACES 
PYEGeEW LSE CONSTANT 


D 
ED 
ATA 


65 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 -~46.0840542053 

2 -40.0840542053 

3 ~3-0000451598 

4 — FeSO LOC92 7 
LEROES REAL 

1 ~19. 9584019372 
GAIN CONSTANT = 160.7285003662 


SYoTEM STATE WARITABLES (PHASE FORM) 


A VECTOR 
if S°1T00.1743/ 
2 19844. 88825 
3 2525 eUSTo52 
4 91.16786366 
B VECTOR 
4 166.7285004 
C VECTOR 
1 19.-95840194 
2 1.-0069C00000 


INITIAL CONDITION (G) VECTOR 


Ik 62.00114277 
2 DOese 549513 
3 160935.4623 
4 1176687.962 


66 


IMAGINARY 
-~4.7409985601 
4-1409985601 
0.0 
0.0 


IMAGINARY 
0.0 





EXAMPLE 9 


Srorem= Vowme ITOENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 -~1.0000000000 

2 -5-9000000000 
ZEROES REAL 
GAIN CONSTANT = 10.0000GG000C 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 52900000060 

2 §- 0000000GCQ 
B WECTOR 

2 16.0099)000CC 
C VECTOR 

l 1.900909C3000 


PNETITAL STATE VECTOR 
1 1.00000C00G0 
- 
2 0.0 


67 


IMAGINARY 
0.0 
C0 


IMAGINARY 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


ECEGS REAL 

i ~5.0110029817 

2 —-28 74. 402942 1/054 

3 C0.00945131 73 

4 C.0094513173 

5 =—() 991999 181363 
BERUES REAL 

1 0.0096691675 

2 0.0096691675 
GAIN CONSTANT = 28777.20156250000 


eect EM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 84255. 60063 
2 1068286134 
3 30926. 28200 
4 17234. 39503 
5 2880.395021 
B VECTOR 
5 28777-01563 
C VECTOR . 
1 54855754701 
2 -0. 1933833505D-01 
3 1.000000000 


68 


IMAGINARY 
0.9 
0.0 
-~224186034399 
2¢4186034399 
0.0 


IMAGINARY 
—2¢ 4198473522 
264198473522 


a. a. a ce CG 





IDEN 


TIFICATION OF UNKNOWN SYSTEM 


INITIAL CONDITION (G) VECTOR 


1 


mm ££ W NY 


1-450 317645 
28764228032 
17252. 85740 
16462.65716 
101232.6568 


PROGRAM PARAMETERS 
RECORD LEGNTH 


= 12.000G0O SEC 


SAMPLING PERIOD = 2e000CO0 MSC 
EQUATIONS FORMED = wae 

PRECISION OF DATA = 5 DEGRMAL PLACES 
INPUT FUNCTION = PIECEWISE CONSTANT 


69 


IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 =e 9999961462 

2 ~5.6000869786 
ZEROES REAL 
GAIN CONSTANT = 10.0001306534 


Seovcem STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 5.-000067709 
2 62000083125 © 


2 10.€C0013065 
C VECIOR 
1 1.000000000 


INITIAL CONDITION (G) VECTOR 
i 1.000025784 
2 6-000036891 


70 


IMAGINARY 
0.0 
0.0 


IMAGINARY 





AMELIE 10 


eroteme 1G BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL 

1 -5.000690C000 

2 -5.0000000000 

3 -100.00900C60000 
ZEROES REAL 

1 ~20.0000000000 
GAIN CONSTANT = 109.0CO00CO0GG00 


Srol GM STATE VARIABLES (PHASE FORM) 


A VECTOR 
Ji 5000.00000C0909 
2 1050. 9000000000 
3 110.009000000C 
B VECTOR 
3 10020099000000 
C VECTOR 
1 20.090000C0000 
2 1-9000000000, ~ 


INITIAL STATE VECTOR 


1 2eCCOCO000000 
Z 12.90000000000 
3 Od 


71 


IMAGINARY 
-5.-0000C909000 J 
52-900000C000 J 


0.0 J 
IMAGINARY 
C.0 J 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES 
l 


wm £ WwW NH 


ZEROES 
1 
(a 


GAIN CONSTANT 


SYSTEM STATE VARIABLES 


A VECTOR 
1 


mn & W NH 


B VECTOR 
5 


Ce VECTOR 
1 
2 
3 


REAL 
3533. 8015984749 
—-252.9614103219 
~100.0395044818 
=“eo7Je911 lol 
=4. 999608 (161 


REAL 
—-19.5630088534 
~26.9189879242 


=~348534. 3 7590009000 


~458893162.9 
=lisgls 225.62 
—-137174844. 37 
~476693.1811 
35916600198 


~348534.3750 


526.26163991 
46248199678 
1. 9000000000 


7 


(PHASE FORM) 


IMAGINARY 
0.0 
0.0 
0.0 

-5-¢0000893933 

5e 0000893933 


IMAGINARY 
C.0 
0.0 


IDENTIFICATION OF UNKNOWN SYSTEM 


INITIAL CONDITION (G) VECTOR 


1 40.97649863 
2 =o Oo 4a o2 
3 =19I623611-4- 16 
4 = SS o7 Stor (24 
5 -3138782043. 


PROGRAM PARAMETERS 
RECORD LEGNTH 
SAMPLING PERIOD 
EQUATIONS FORMED 
PRECISION OF DATA 
INPUT FUNCTION 


C.6000C SEC 
O.10CCG MSC 

120 

> DECIMAL PLACES 
PIECEWISE CONSTANT 


13 





IDENTIFICATION OF UNKNOWN SYSTEM 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 S997 9 OO O96 2 0.0 J 
2 ~4.°9999968528 ~520000354476 J 
3 = aes 2 OD. 2 6 520000354476 J 
ZEIGE S REAL IMAGINARY 
1 ~19.9960841378 0.0 J 
GAIN CONSTANT = 100.0136718750 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 4999.8 30647 
2 1649.959363 
3 1 96 SIDS Of 
B VECTOR 
3 100.0136719 
CeVECTOR 
1 1999608414 
@ 1. 0600990090 


INITIAL CONDITION (G) VECTOR 


1 49.99999105 
Ze 45292627485 
S 34198. 6335C 


74 





IV. CONCLUSIONS 


The method of multiple integrations is a Praceil cal 
and flexible method for identifying lumped parameter, 
linear, time invariant systems. Complete and accurate 
identifications can be made on the basis of arbitrary in- 
put-output records taken over a short time interval and 
pecuvate to only thmee or four significant figures. The 
computational requirements of the method are not exces- 
Sive. The procedure can be implemented by relying en- 
tirely on subroutines which are available in most computer 
Center libraries. 

At present the technique is limited to comparatively 
low order systems. When the input-output records are good 
to three or four significant digits the method will be 
capable of identifying systems up to about fifth order. 
This limitation is due primarily to the algorithm used he: 
solve the overdetermined set of linear equations. As 
better algorithms become available it will be possible to 
identify higher order systems. 

The accuracy of an identification is usually compar- 
able to the number of significant digits in the Lnpwue- 
output data. Accuracy depends to a lesser extent on the 
Sampling rate, the nature of the input function dri ime 
the system, and the order of the system. 

There are several areas where additional research 


might prove fruitful. Since the input-output records must 


iis 





be in sampled data form for the computer it might be prof- 
meable to neformulate the identificationsprocedure from a 
sampled system standpoint. Standard techniques could be 
used to convert the sampled system representation obtained 
from the identification to a continuous system representa- 
tion. The ability to identify sampled systems would be a 
worthwhile extension of this method. 

The method of multiple integrations may prove very 
useful as a tool for approximating high order systems with — 
low order systems. Research could be done to determine 
the quality of the approximations obtained using this 


method. 
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AME 


Mm OAMNM 


MAM AMNAOY 


OOO 


OOO 


NAD OAM 


OOO 


10 


15 


20 


Zo 


30 


35 


40 


45 


50 


NZ1=NZ+1 

N=NP+NZ41 

NI=N+NP 

SET CUMULATIVE INTEGRAL VALUES TO ZERO 
DO 10 1=2,NP2 

RO(1)=0.00 

CO(1)=0.0C 

READ HN INITIAL DATA POINT (1)R>C) 
CALL SYSTEM(TO,RO(1),CO(1)) 


TOFF=TO 
K=K+1 


READ IN NEW DATA POINT (T,RsC) 
CALL SYSTEM(TN(1),RN(1),CN(1)) 
UPCATE MULTIPLE INTEGRATIONS 


OT2=(TN(1)~-TO) *0.5 

DO 20 INT=1,NP1 
RNCINT+1I=CROCINTI +RNCINT) )*¥OT2+ROC INT +1 ) 
CNCINT41)=(COCINT) +CNCINT) ) *¥OT24+COCINT+1) 
FORM A LINEAR EQUATION 

eee NEs tier IPTS)*IPTS) GO TO 35 
M=M+) 

B(M)=CN(2) 

TN(2)=(TNC1)-TOFF) 

DO 25 T=1;NP 

TA=(NP-T)*MEQS+M 

Pe= (+ I-] )*¥MEQS +h 

WNC 2) =FNC14+1)* TNC 2) /FLOATCI+1) 
AC IA)=-CN(1+2) , 
ACIC)=TN(IT+1) 

DO 30 1=1,NZ1 

TA=(NP+I~-1)*MEQS+M 

IRN=NP+3-I 

A(TA)=RNCIRN) 

Rieeer OCD VALUES 

TO=TN(1) 

DO 40 I=1,NP2 

RO(I)=RN(T) 

COC IJ=CN(T} 

Te (MetTeMEQS) GO TO 15 


SOLVE FOR PARAMETERS BY METHOD OF LEAST SQUARES 


CALL DLLSQ(A,ByMyNI 91s XeTPIVyEPSs TER sAUX) 
CRECUEATE POLES 


AP(1)=1.00 

0 45 1T=1,NP 

=NP+2-I 

(lx (1) 

LL RTPLSB(NPyAPy PRA, PIA;CONV, IERPZ) 
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GO TO 60 
NZ;AZ;ZRA,;ZIA;,;CONV,IERPZ) 


omy ww 


OUTPUT 
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EXO oO. Ni ht <I : < 
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SUBROUTINE MINV 


PURPGSE 
INVERT A DOUBLE PRECISION MATRIX BY THE GAUSS- 
JORDAN METHOD. 
USAGE 
CALL MINV(CASN,D,LyM) 
DESCRIPTION OF PARAMETERS 
A - INPUT MATRIX, DESTROYED IN COMPUTATION 
AND REPLACED BY RESULTANT INVERSE. 
N - @EDER OF MATRIX A 
D -= "RESULTANT DETERMINANT 
Ie = weRK VECTOR OF LEGHTR ON 
M - WORK VECTOR OF LEGNTH N 
TION SUBPROGRAMS REQUIRED 


SUBROUTINES AND FUNC 
NONE 

REMARKS 
NONE 


. 
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SUBROUTINE MINVCAN1O,L,M) 
DIMENSTON A(1),L¢1) M01) 
DGUBLE PRECISION A,0;BIGA;,HOLD 


SEARCH FOR LARGEST ELEMENT 


D=1.0 

NK=—N 

00 80 K=1,N 

NK=NK+N 

L(K)=K 

M(K)=K 

KK=NK+K 

BIGA=A( KK) 

DO 20 J=K,N 

TZ=N*(J-1} 

DO 2G I=K,N 

[Jje17+1 

IF (DABS(BIGA)-DABS(A(IJ))) 15,20,20 

BIGA=A( IJ) 

L(K)=I 

M(K)=J 

CONTINUE 
INTERCHANGE ROWS 

J=L(K) 

IF(J-K) 35,35,25 

KI=K-N 

DO 30 T=1,N 

KI=KI4+N 

HOLD=-A(KI) 

JT=HKI~K+3 

AC(KI)=AC( JIT) 

A( JI) =HOLD 
INTERCHANGE COLUMNS 

I=M(K) 

IF(I-K) 45,45, 38 

JP=N*(I-1) 

DO 4G J=1,N 

JK=ENK4F J 

JI=HJP+H+I 
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SUBROUTINE PROD 


PURPOSE 

TO CGMPUTE THE PRODUCT OF TWO MATRICES 
USAGE 

CALL PROD(A,B;CyNyMyL) 
DESCRIPTION OF PARAMETERS 

A —- FIRST INPUT MATRIX 

B ~ SECOND INPUT MATRIX 

6 - PRODUCT OF A AND B 

N ~ NUMBER OF ROWS IN A AND C 

M ~ NUMBER OF COLUMNS IN A AND ROWS IN B 

L ~ NUMBER OF COLUMNS IN B AND C 
REMARKS 

NONE 


Sc, ae AND FUNCTION SUBPROGRAMS REQUIRED 


E PROD(A,ByCyNyMyL) 
19) 9B(09;9) +C(949) 
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SUBROUTINE SUMM 


PURPOSE 
ADD TWO MATRICES 
USAGE 
CALL SUMM(A,B,CyM,N) 
DESCRIPTION OF PARAMETERS 
A — Naame OF FIRST INPUT MATRIX 
B ~ NAME OF SECOND INPUT MATRIX 
C ="Nemt OF OUTPUT MATRIX 
M ~ NUMBER OF ROWS IN AyBy AND C 
N - NUMBER CF COLUMNS IN AyBy AND C 
REMARKS 
NONE 


ee AND FUNCTION SUBPROGRAMS REQUIRED 
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SUBROUTINE SUMM(A,B,CeM,N) 


REAL*8 A(9,9),B(9;9),C(9,9) 
DO 1 T=1,M 

DO 1 J=1,N 
C(I,JI=AC1I,J)+B801,J) 

RETURN 

END 
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SUBROUTINE DIFF 


PURPOSE 
SUBTRACT ONE MATRIX FROM ANOTHER 


USAGE 
CALL DIFF(A,ByCyMyN) 


DESCRIPTIGN GE. PARAMETERS 

A = Fis ere MATRIX 

B - SECOND ENPUT MATRIX 

C ~ GUTPUT MATRIX EQUALS A — B 

M - NUMBER OF ROWS IN A>B, AND C 

N ~ NUMBER OF COLUMNS IN A,yBy AND C 
REMARKS 

NONE 


oS ee AND FUNCTION SUBPROGRAMS REQUIRED 


SUBROUTINE DIFF(A;B,CeM;N) 
REAL¥8 A(9,9),B(9;, Sap 9) 
DO 1 Y=1,M 

DO 1 J=1,N 

C(I,JJ=HA(I, JI-B(I, J) 

RETURN 

END 
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SUBROUTINE EXPAND 


PUREOSE 
WerCOMPUTE THE REAL COEFPFICUENTS GPaAnei a 
DEGREE POLYNOMIAL GIVEN N COMPLEX ROOTS 
USAGE 
CALL EXPAND(N;R,A) 
DESCRIPTION OF PARAMETERS 
N DEGREE OF POLYNOMIAL 
R -~(VECTOR CF COMPLEX Tears 
A —-V COEFFICIENT V eerie 
REMARKS 
NONE 


SS aoa AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 


ASE aN EXPAND(N;R;A) 
*8 A(10) 
EX*¥16 R(9),Q(9) 
Pele Liye s 
e O00 
REAL(R(1)) 
2 OF 
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